The role of long waves in the stability of the plane wake 
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This work is directed towards investigating the fate of three-dimensional long perturbation waves 
in a plane incompressible wake. The analysis is posed as an initial-value problem in space. More 
specifically, input is made at an initial location in the downstream direction and then tracing the 
resulting behavior further downstream subject to the restriction of finite kinetic energy. This presen- 
tation follows the outline given by Criminale and Drazin [Stud, in Applied Math. 83, 123 (1990)] 
that describes the system in terms of perturbation vorticity and velocity. The analysis is based 
on large scale waves and expansions using multi scales and multi times for the partial differential 
equations. The multiscaling is based on an approach where the small parameter is linked to the 
perturbation property independently from the flow control parameter. Solutions of the perturbative 
equations are determined numerically after the introduction of a regular perturbation scheme ana- 
lytically deduced up to the second order. Numerically, the complete linear system is also integrated. 
Since the results relevant to the complete problem are in very good agreement with the results of 
the first order analysis, the numerical solution at the second order was deemed not necessary. The 
use for an arbitrary initial-value problem will be shown to contain a wealth of information for the 
different transient behaviors associated to the symmetry, angle of obliquity and spatial decay of 
the long waves. The amplification factor of transversal perturbations never presents the trend - a 
growth followed by a long damping - usually seen in waves with wavenumber of order one or less. 
Asymptotical instability is always observed. 



PACS numbers: 47.15.Fc, 47.15.Tr, 47.11. St, 47.20.Ft 



I. INTRODUCTION 

The traditional investigation of stability of shear flows 
is cast as a linear initial-value perturbation problem. In 
principle, save for the additional complexity of necessitat- 
ing three space dimensions as well as time, this is done by 
means of a Laplace transform in time. Once the bound- 
ary conditions have been satisfied, the stability or non 
stability is found. Further, depending upon the mean 
shear flow that is being investigated, the causes are deter- 
mined. No attention is given to any specific input or the 
effect of various physics. Moreover, little attention was 
given to early period dynamics, see for example Grosch 
and Salwen Salwen and Grosch Q. These authors 
showed that there can be early time growth of a pertur- 
bation even if there is damping for long time. In short, 
a branch cut can be present as well as any pole when 
inverting the Laplace transform. Next, from the labo- 
ratory, interest turned to spatial growth or decay after 
an input at an initial location rather than the temporal 
behavior. This construction creates new difficulties but 
they are not insurmountable [|, I|. Stih, just as in the 
temporal problem, no specific initial input has been ex- 
amined. Regardless of the framework, it has been known 
since the first results of stability theory, that the value 
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of the wavenumber that comes into the analytical frame- 
work (due to Fourier decomposition in the variables in 
the plane that is perpendicular to the mean flow) is small 
in the regions where there is instability. In short, long 
waves. Such a result provides a sound means for the anal- 
ysis and examination of a specific initial input. This is 
true whether posed as a temporal or spatial initial-value 
problem. It further provides a means to investigate in- 
teraction, the early period and a way to suppress any 
growth at the early period or location. In 1962 a study 
about the instability to long waves of unbounded par- 
allel inviscid flow was given by Drazin and Howard [Bj. 
Using the normal mode analysis, they found that there 
is a finite number of different modes unstable to long 
waves, essentially one for each relative maximum or min- 
imum of the velocity profile. Healey considered long 
waves for investigating spatial instability of the rotating- 
disc boundary layer, and by means of an analytic theory 
in the inviscid long wave limit, he obtained an explicit 
expression for the growth rate in terms of basic flow pa- 
rameters. 

Large or long waves have now been used in full non- 
linear simulations. For example, Ryzkov and Shcvtsova 
focused on convective instability in multicomponent 
fluids, showing by means of both linear stability analysis 
and nonlinear numerical calculations that the instability 
is caused by the interplay between the basic flow and the 
concentration waves which have a long scale in a vertical 
direction. And Barros and Choi la considered the inhibi- 
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tion of the shear instabihty that can be induced by large 
amplitude internal solitary waves travelling in a two-layer 
flow with a top free surface. For large eddy simulations 
in turbulence see @ or [loj among many others. 

The analysis in the present work is based on large scale 
waves and expansions using multi scales and multi times 
for the partial differential equations. The multiscale is 
based on an approach where the small parameter is linked 
to the perturbation property independently from the flow 
control parameter. In fact, the perturbation scheme is 
based on the introduction of a small parameter which is 
the wavenumber k in the limit fc — >■ 0, and is analytically 
determined to the second order. The perturbative equa- 
tions used follow the formulation given by Criminale and 
Drazin fllj that describes the system in terms of pertur- 
bation vorticity and velocity. Numerically, the complete 
linear system has also been integrated for the non paral- 
lel base flow. What results is an extension of a previous 
work based on a locally near parallel assumption (T^j. 
The formulation of the linear perturbation initial-value 
problem is presented in Section 2. Results are in Section 
3. Conclusions follow in Section 4. 



II. FORMULATION 

By exciting the plane wake flow 
{\J=(U{x,y; Re)^V{x,y; Re))) with small arbitrary 
three-dimensional perturbations, the continuity and 
Navicr-Stokcs equations for the perturbed system 
linearized with respect to small oscillations are given by 
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where (u{x,y,z,t), v{x,y,z,t), w{x,y,z,t)) and 
p{x^y,z,t) are the components of the perturbation 
velocity and pressure, respectively. 

The independent spatial variables z and y are defined 
from — cx) to +oo, x from to +oo. All physical quantities 
are normalized with respect to the free stream velocity 
[//, the body scale D and the density p. The Reynolds 
number is defined as Re = UfD/v, where v is the kine- 
matic viscosity. 

The two dimensional wake is a thin free flow that can 
be schematized as shown in FigHJ Leaving aside the 
near field, that is highly non parallel since it hosts the 
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FIG. 1: Base flow sketch. The base flow has been chosen in 
order to be an acceptable representation of the intermediate- 
far fleld. To this aim we build a homogeneous fleld in x, z 
by using the information associated to a section, xo, placed 
in the intermediate region, x £ [5,30]. In the sketch the 
longitudinal and transversal proflles at Re = 100 are frozen at 
xq = 10 (note that the transversal velocity V is multiplied by 
a factor 10). The base flow {U{y; xq, Re), V{y; xq, Re)) is thus 
a slightly non parallel flow homogeneous in x,z, which makes 
it possible to Laplace transform the perturbative equations in 
X and to Fourier transform them in z. 



two symmetric counter circulating vortices that consti- 
tute the separation region, the intermediate and long 
term wake is a near parallel flow. The wake slowly be- 
comes thicker according to a law which, at first order, 

X ~ 

scales as ( — ) . As representation of this steady sub- 

critical flow we consider the asymptotic expansion so- 
lution in inverse powers of x obtained in [13|. In par- 
ticular, we consider the intermediate far field well rep- 
resented by a section placed near x = xq = 10 and 
build the basic flow by freezing it at this longitudinal sta- 
tion. In so doing, the basic flow is parameterized through 
the downstream station xq and the Reynolds number Re 
{lJ={U{y; xo,Re), V{y; xq, -Re))). It is thus homogeneous 
in X and z. As such, the long waves that are the main 
subject of this study are valid. 

The explicit expressions of the base flow longitudinal 
and transversal components are: 



U{y]Xo,Re) = + (j)iXf^ 

I J, -1 1 J, -3/2 
— 1/2 

V{y;xo,Re) = Xo + XiXq 



I -1 I -3/2 

+X2X0 +X3X0 



(5) 



(6) 



The coefficients (pi = (j)i{y;xo, Re) and Xi = 
Xiiy, Xq, Re) of this expansion up to i = 3 arc given 
in Appendix A. Fig. 2 displays the intermediate wake 
profile for which there exists a comparison based on 
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FIG. 2: Example of velocity profiles in the intermediate wake 
at the downstream station x — xo = 20. (a) Longitudinal 
velocity U, (b) transversal velocity V. Continuous curves: 
analytical solutions {Re = 20, 60 and 100) by Tordella and 
Belan (2003) [13 [.triangles: numerical results {Re = 34) by 
Berrone (2001) 14[]. circles: laboratory data {Re = 34) by 
Kovasznay (1948)^. 



laboratory and numerical simulation results, 

fil,[il,fii, M. 



see also 



It should be noted that when using such a kind of rep- 
resentation the base flow nonparallelism is considered and 
allows for the effect of the lateral entrainment to be ob- 
tained [l6j . However, in this study a fixed location xo of 
the intermediate wake is considered since in this region 
absolute instability pockets have been found by recent 
modal analyses [l^, [13 • The term intermediate is used 
in the general sense as that given by Barenblatt [l^: 
intermediate asymptotics are self-similar or near-similar 
solutions of general problems valid for times and dis- 
tances from boundaries, large enough for the influence 
of the fine details of the initial or boundary conditions to 
be insignificant, but small enough that the system is far 
from the ultimate equilibrium state...'. The distance be- 
yond which the intermediate region is assumed to begin 
varies from eight to four diameters D for Re e [20, 40] 

m. 
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where {ujx,ujy,ujz) is the perturbation vorticity field, 
dV dU 

^lz = (-;: K~)\x=xn is the mean vorticity in the span- 
ax dy 

wise direction, and the coefficients a, 6, c, d, e are the 
spatial derivatives of the base flow vorticity and velocity 
at Xo, namely: 



dn. 



dx 
dU 
dx 



dx^ 



c(y) 



dxdy 



X — Xq 



dV 
dx 



X — Xq 



By introducing the quantity F, that is defined by 

f (9) 



we obtain three coupled equations ([7]), ([8]) and ([9]). Equa- 
tions ([7]) and (|S]) are the Orr-Sommerfeld and Squire 
equations respectively, from the classical linear stability 
analysis for three-dimensional disturbances. From kine- 
matics, the relation 
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(10) 



physically links the perturbation vorticity components in 
the X and z directions {oJx and u)z respectively) and the 
perturbed velocity field. By combining equations ([7]) and 
dH) then 
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which, together with (|5]) and (jH]), fully describes the 
perturbed system. Since seven unknown quantities 
(m, V, w, uixi^Oy^oOz, r) are involved in the above equations 
dS]), ([U and ((TT|) . the perturbation vorticity definition 
and the continuity equation 



By combining the momentum equations to (j?]) to 
eliminate the pressure, the resulting governing equations 
become 



LO — V X U, 
V - M = 0, 
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(13) 
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FIG. 3: Perturbation geometry scheme. 



FIG. 4: The wave spatial evolution in the x direction for k : 
Or = 0.05, Qi = 0.1,0.01,0.001. 



link the pcrturbative system of equations ([5]), © and 

m- 

For every dependent variable, we perform a combined 
spatial Laplace-Fourier decomposition in the x and z di- 
rections defined by 



+ 00 /' + 00 



g{x,y,z,t)e~ 



^dxdz 



where 5 is the general dependent variable, a, the lon- 
gitudinal wavenumber, is complex {a = Ur + iai) and 
7, the transversal wavenumber, is real. By adopting 
the velocity- vorticity formulation, [ll|, Q, the govern- 
ing equations (|5]), ^ and ([TT]) can now be written as 



dt 
'dt 

dtUy 

'df 



(fc^ -af + 2ikcos{(j>)ai)v = t (14) 



— = GT + Hv + Kcjy (15) 



Luiy + Mv (16) 



where cf) = tan^^{'j/ar) is the perturbation angle of 
o bliquity with respect to the x-y physical plane, k = 
y^a^~\-j^ is the polar wavenumber, at is the imaginary 
part of the complex longitudinal wavenumber, ujy is the 
transversal component of the perturbation vorticity, and 
F is the vorticity component in the oblique direction 
which is defined as F ~ i{cti^z — I'^x)- In Figure [3] 
the three-dimensional pcrturbative geometry scheme is 
shown. 

In order to have a finite perturbation kinetic energy, Ui 
can only assume non-negative values. In so doing, we 
allow for pcrturbative waves that can spatially decay 
(ckj > 0) or remain constant in amplitude (a,; = 0). In 
the following, Ui is called spatial damping rate. It should 



be pointed out that the present analysis is not the stan- 
dard eigenvalue problem where poles result. Here, in fact, 
the spatial damping rate ai is a parameter and, as such, 
should be simply imposed. The magnitude of the spatial 
damping rate can vary in order to describe a physically 
meaningful damping of the pcrturbative wave in the x 
direction (disturbances immediately damped to zero are 
not allowed). According to this and since long waves 
[k ^ 10^^, 10^^) are considered, ai is non-negative and 
at maximum ~ 10""'^, see Fig. 01 Symbols G, H, K, L 
and M represent ordinary differential operators, written 
in the form G = G{y;xQ,k,(j),ai, Re), and similarly for 
H, K, L and M, since they arc functions of y and are 
parameterized through the fixed longitudinal station xq, 
the polar wavenumber k, the angle of obliquity 0, the 
spatial damping rate ai and the Reynolds number Re. 
All these operators are explicitly given in Appendix B. 

Equations (jTS]) and require proper initial 

and boundary conditions in order to be solved. As far 
as the boundary conditions are concerned, among all so- 
lutions, those whose perturbation velocity field vanishes 
in the free stream are sought. The initial conditions are 
necessary for F and ujy. As far as the initial conditions 
for F are concerned, according to equation (6), they can 
be shaped in terms of set of functions in the Hilbert 
space via the variable v, which is here represented by the 
trigonometric system 



v{0,y)^e ^ cos(y), v{0,y) = e ^ sm{y), 

for the symmetric and the asymmetric perturbations, re- 
spectively. This trigonometrical system is a Schauder 
basis in each space Lp[0, 1], for 1 < p < 00. The transver- 
sal vorticity LUy is instead taken initially equal to zero in 
the y domain, in order to directly observe the net con- 
tribution of three-dimensionality on its temporal evolu- 
tion. However, it can be demonstrated that the eventual 
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introduction of an initial transversal vorticity docs not 
substantially affect the perturbation temporal evolution. 

In the stability analysis of spatially developing flows, 
different scales can be determined. Usually, long and slow 
scales, related to the slow base flow evolution, as well as 
short and fast scales, linked to the disturbance dynam- 
ics, can be defined. However, it should be noticed that in 
some flow configurations, long waves can be destabiliz- 
ing. Examples of this behavior are the two-dimensional 
Blasius boundary layer, the three-dimensional cross-flow 
boundary layer, as well as the free shear flows. In such in- 
stances, the perturbation wavenumber is less than 0{1) 
where instability occurs. Thus, a regular perturbation 
scheme can be adopted, defining the polar wavenumber 
k as the small parameter [l^], [2l|. It should be noted 
that by using such a long-wave expansion, the x-scale 
length of the perturbations is comparable to the x-scale 
length of the base flow. Indeed, wc only consider the in- 
termediate and far wake sections, where the flow slowly 
evolves in the longitudinal direction. Thus the near wake 
is not taken into account. In synthesis, the scale of the 
intermediate wake is of the order xq ^ 10^, and the scale 
of the long pcrturbative waves is A = 2TT/k ^ 10^, 10^. 

Two spatial scales, a short one, y, and a long one, 
Y ~ ky, are defined. For the temporal dynamics, three 
temporal scales, the fast one, t, and the slow ones, t = kt 
and T — k^t, can be identified. The perturbation quan- 
tities {v, r, tby) are now function of y, Y, t, r, T, expressed 
as r = r(y, y, t, t, T; fc, (/), tti), and similarly for v and ujy, 
and can be expanded as 
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f 



k'^V2 



vo + kvi 
fo + fcfi + k^V2 + 

LJyQ + k(Jj. 



(17) 



Initial conditions at order 0(1) are defined as in the full 
linear problem, while at higher orders (0(fc), 0{k^), . . .) 
are equal to zero. Boundary conditions remain as stated 
in the full linear problem. Substituting relations (jl7|) in 
the full linear system (|14p to (|T6l) . the following ordered 
hierarchy of equations, expressed up to 0{k), result and 
are: 

Order 0(1) 



duj. 



dt 



- HqVo 

Lo(jJyO 



(18) 

(19) 
(20) 



dt 



dydY 
Gof 1 - i/owi = 
5fo 



2icos{(j))aiVo -j- Ti (21) 



dr 

LoLOyl = 

~ dr 



GiTo + Hivo + Kiiuyo (22) 



(23) 



Order 0(k) 



Operators Gq = Go{y;xQ,(l),ai, Re) as well as and 
Lq are functions of the short scale y only. Operators 
Gi = Gi{y, Y; xq, (j), a^. Re) as well as Hi, Ki, Li and Mi 
are function of both the short scale y as well as the long 
scale Y. These operators are explicitly given in Appendix 
B. 

A comment concerning the role of ai is needed. Equa- 
tions (10) to (15) above are obtained for the case where 
ai > 0. It can be observed, see the Appendix B, that 
if cti = 0, the 0(1) operators Ho,Lq (but also the 
0{k) operators Hi, Li, Mi, Ni) are singular. It is pos- 
sible to verify that if a; = the disturbances initially 
imposed remain constant as time passes and reach, in 
the end, an asymptotic condition of marginal stability. 
This fact is deduced by considering equation p8)) . For 
ai = 0, the homogeneous solution assumes the expression 
voh = ci + C2y. Since the perturbation velocity field has 
to vanish in the free stream, ci = and C2 = 0. Thus 
vq, and therefore Fq, vanish as well. This means that, 
in equation (fT9| . = 0, so that there is no temporal 
evolution for Fq. Since the transversal vorticity ujyo is 
initially taken as zero, then in equation ^U^, = 0, 
and thus, also for the transversal vorticity, there is no 
temporal evolution. The complete problem is defined for 
ai = and, for this value of ai, it does not necessarily 
show a condition of marginal stability (e.g. sec Fig. 6 
in the following). However, it is possible to see that the 
multiscaling limit for — )■ well approximates the com- 
plete problem: cf. Figures 5 and 6. When ai = 0, the 
multiscaling presents a discontinuity, since it has a right 
limit different from the value shown at = 0. 

Order 0(1) is the more important approximation of the 
pcrturbative analysis and its formal expression is simpli- 
fied with respect to the complete problem. Numerically, 
the complete linear system was also integrated. Since the 
results relevant to the complete problem arc in very good 
agreement with the results of the first order analysis, in 
the present work, attention is focused on the resolution 
of the multiscaling at order 0(1). It should be noted that 
at this order only the short spatial scale y is of relevance. 

In the following, a summary of the most signifi- 
cant transient behavior and asymptotic fate of three- 
dimensional perturbations is presented to highlight the 
agreement between solutions of multiscaling at order 
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0(1) and full linear problem. Results will be principally 
focused on parameters such as the spatial damping rate, 
the polar wavenumber value (to check the validity of the 
approximation), the angle of obliquity and the symmetry 
of the three-dimensional disturbance. 

To measure the transient growth the concepts of ki- 
netic energy density e(t; fc, (j), a^) 



e(t;fc,(/>,a,) = i-^ / ""'Hul^ + \v\^ + \w\^)dy 
2 2yd J^y, 



1 1 



1 



X 



(24) 



2 2yd \ P + 2ikcos{4))ai ~ af\ 
r+Vd 

/ + \k'^ + 2ikcos{(t>)a^^ a^,\\v\'^ + \u;y\^)dy, 

J-ya 9y 

and normalized amplification factor G(t; fc, 0, aA 
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G{t; k, (j), a,) = -J- — — ^. 25) 

e[t = 0; k, (j), ttj) 

are introduced for both multiscale 0(1) quantities 
(vo,TQ,iliyo) and full problem quantities {v,T,iliy). 

In (16), the limits ±yd define the spatial extension of 
the numerical domain. The value yd is defined so that the 
numerical solutions are insensitive to further extensions 
of the computational domain size. Here, in the limit of 
long waves, the size of the spatial domain 2yd assumes 
values in the range between 30 and 100 external flow 
scale D. The total kinetic energy can be obtained by 
integrating the energy density over all k and (f>. 

To evaluate the asymptotic behavior we introduce the 
temporal growth rate r, defined as 

I, A \ ^og|e(^;a,7)| 
r{t;k,(j),a.,) = — , t>0. (26) 

The temporal growth rate r is not defined for t = 0. This 
quantity has, in fact, a precise physical meaning asymp- 
totically in time. Moreover, for both multiscale and the 
full problem solutions, the angular frequency (pulsation) 
uj of the perturbation can be introduced by defining a 
local, in space and time, time phase (f of the complex 
wave 

viy, t; a, 7, Re) = At{y; a, 7, i?e)e^'^(*) (27) 

and then computing the time derivative of the phase per- 
turbation ip 

, , , , X dip{y,t;k,(j>,ai) 
uj{y, t; k, cj), a,) = — . (28) 

Since Lp is defined as the phase variation in time of the 
perturbative wave, it is reasonable to expect constant 
values of frequency, once the asymptotic state is reached. 



■ 10 20 t 30 40 50 

FIG. 5: Effects of the spatial damping rate a^. (a) The 
amplification factor G and (b) the temporal growth rate r 
as function of time. Comparison between multiscale 0(1) 
(thick curves) and full problem (thin curves). Re — 50, 
k = 0.03, (fi = n /A, xq — 12, asymmetric initial condition, 
ai = 0.04, 0.4. 



III. RESULTS 

Computations to evaluate the long time asymptotics 
are made by integrating the equations forward in time 
beyond the transient until the temporal growth rate r, 
defined in relatio n ([261) . asymptotes to a constant value 
{dr/dt < e) m, [lip 

Fig. [S] presents an interesting phenomenon that is ob- 
served for general long perturbations (either transversal, 
or oblique, or longitudinal) by changing the value of the 
spatial damping ai. For instance, in the case shown in 
this figure, which is relevant to a long oblique asymmetric 
wave, the variation of the order of magnitude of ai from 
0.04 to 0.4 highly enhances the amplification in time, 
with a temporal growth rate that becomes nearly three 
times larger. This means that perturbations that are 
spatially confined are more amplified in time. It can also 
be noted that the agreement between multiscale 0(1) 
(thick curves) and full problem (thin curves) remains very 
good when changing the order of magnitude of the spatial 
damping. 

The influence of the perturbation symmetry on the 
early time behavior is shown in Fig. [5] (a logarithmic scale 
is used on the ordinate of part (a) of the figure). It can 
be noted that the symmetric initial condition leads - in 
the transient behavior - to a faster temporal growth than 
the asymmetric one, although both configurations are ap- 
proaching the same asymptotic unstable state. Indeed, 
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FIG. 6: Transversal long waves. Effect of the symmetry of 
the perturbation. Comparison between multiscale 0(1) (thick 
curves) and full problem (thin curves), (a) The amplification 
factor G and (b) the temporal growth rate r as function of 
time. Re = 100, k = 0.02, <j> = tt/2, xq = 10, = 0.08, 
symmetric and asymmetric initial conditions. 



FIG. 8: Transversal long waves. Effect of the spatial decay, 
Oi. Comparison between multiscale 0(1) (thin curves) in the 
limit for — ^ 0, and full problem (thick curves) with Oi = 0. 
(a) The amplification factor G and (b) the temporal growth 
rate r as function of time. Re = 50, (j) = 7r/2, xq ~ 10, 
asymmetric initial condition, k = 0.04, Oi = 0.005,0.01,0.05 
(multiscale 0(1)), Qi — (full problem). 
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FIG. 7: Effects of the polar wavenumber k used as small 
parameter, (a) The amplification factor G and (b) the tem- 
poral growth rate r as function of time. Comparison between 
multiscale 0(1) (thick curves) and full problem (thin curves). 
Re = 100, (j> = 0, xo = 27, Oi = 0.2, symmetric initial condi- 
tion, k = 0.1,0.01,0.001. 



FIG. 9: Temporal asymptotic values of : (a) the tempo- 
ral growth rate and (b) the angular frequency. Comparison 
between multiscale 0(1) (squares: symmetric inputs, dots: 
asymmetric inputs) and full linear problem (circles: symmet- 
ric inputs, triangles: asymmetric inputs), k = 0.01, Re = 100, 
= tt/4, xo = 10, dr/dt < e with e ~ 10"*. 
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the transient in the asymmetric case is lasting longer than 
t ~ 10^, where in the symmetric case is lasting t ^ 10^. 
The agreement between multiscale to 0(1) and the full 
problem is very good for both asymmetric and symmet- 
ric conditions. This is true both for the early transient 
as well as the ultimate fate. It should be noted the dis- 
continuous behavior shown at t ^ 100 by the tempo- 
ral growth rate of the asymmetric transversal {(f) = 7r/2) 
wave. 

The effect of differing orders of magnitude for the polar 
wavenumber k is highlighted in Fig. [71 Three orders are 
considered, namely k = 0.1,0.01,0.001. As expected, for 
smaller values of the polar wavenumber the agreement 
between multiscale 0(1) and full problem is improving 
(the multiscale 0(1) solution practically coincides with 
that of the full problem for k = 0.001). It is interesting 
to note in the interval between zero and forty base flow 
time scales, the presence of a temporal oscillation char- 
acterized by a period of about 10 time scales. The tran- 
sient thus presents a further time scale beyond its proper 
global one which, in this case, is lasting 100 units. 

The limit for a zero spatial decay, i.e. — )■ 0, was 
considered in different situations, see Fig. [SJ transients 
of an orthogonal long wave perturbation, and Fig. ^ 
time asymptotics of an oblique long wave perturbation. 
In Fig. [U the thick curves represent the full problem 
solution with Ui = 0, while the thin curves are the mul- 
tiscale 0(1) results with ai values going to zero. The 
right limit of the multiscale 0(1) solution for — )■ is 
finite, and is closely reaching the full problem solution. 
As can be observed, the curves with smaller spatial decay 
rates are approaching the thick curve from above. This 
behavior holds in the early transient and in the asymp- 
totic state. It should be noted that, in this particular 
case (Re — 50, (j) ~ '"'/S, 2^0 — lOi asymmetric initial 
condition, k = 0.04), the complete problem at = 
has a temporal growth rate close to zero, and thus is in 
a near state of marginal stability. But, as previously re- 
marked (cf. Section II), this is not a general behavior for 
the complete problem. However, it is true that the lim- 
iting behaviors for a; — ^ of the multiscaling and of the 
complete problem are very close. And, as the case shown 
in Fig. 6 confirms that, if a difference exists, this will be 
located just at ai = 0. This means that the right limit 
of the multiscaling for a; = correctly approximates the 
complete problem, but this limit value it is not always 
equal to the value at a, = (where the multiscaling 
yields marginal stability, i.e. r = 0, G = 1, w = 0). It 
can be concluded that, in any case, the true limit of the 
complete problem can be obtained by extrapolating the 
multiscaling results. 

It is noted that, in Fig. [8l for asymmetric and transver- 
sal initial conditions with a non vanishing spatial decay, 
a discontinuity in the temporal growth rate can again be 
observed a,t t ^ 100, see also Fig. 3. 

The comparison between the long waves temporal 
asymptotics of the full problem and its multiscaling ver- 
sion is shown in Fig. [HI To consider a situation where the 



multiscaling applies the polar wavenumber k is fixed to 
the value 0.01, while the decay in space at of the longitu- 
dinal wave is in the range [0,0.6]. Multiscale to 0(1) 
results (squares and dots) are in excellent agreement, 
for symmetric and asymmetric initial inputs, with full 
problem data (circles and triangles). Note that the agre- 
mcnt improves for increasing values oi ai. A minimum of 
the perturbation energy (in terms of r) is found around 
ai ~ 0.2 — 0.3 and a similar behavior is shown by the 
angular frequency lj. 

IV. CONCLUSIONS 

Different transient configurations have been observed 
by changing the spatial damping rate, the symmetry of 
the perturbation and its polar wavenumber (the mag- 
nitude, in order to check the method accuracy, and the 
angle of obliquity) . Since the results relevant to the com- 
plete problem are in very good agreement with the results 
of the first order analysis, in the present work attention 
was focused on the resolution of the multiscaling at order 
0(1) only 

Two main results can be noted. First, the perturbation 
symmetry influences the transient. In particular, asym- 
metric transversal perturbations show a different kind of 
transient which includes an initial decay (first few time 
scales) and then a growth that abruptly changes its time 
derivative after about 100 time scales. A sequence of such 
a kind of discontinuities can be envisaged up to where the 
growth rate of the corresponding symmetric perturbation 
is met. Second, the spatial decay substantially affects the 
transient. For example, in the case of asymmetric pertur- 
bations, it is observed that high spatial damping makes 
the initial temporal decay interval shorter and, at the 
same time, greatly increases the temporal growth rate. 

Multiscale data have been compared with full problem 
results in the asymptotic temporal limit. As far as small 
wavenumbers are concerned, the agreement is very good 
for both symmetric and asymmetric initial conditions as 
arbitrarily expressed in terms of elements of the trigono- 
metrical Schauder basis for the space. 

Lastly, it is noted that the amplification factor of 
transversal perturbations never presents the trend - a 
growth followed by a long damping - usually observed in 
waves with wavenumber of order one or less. Asymptot- 
ically unstable configurations in time have always been 
observed here in the limit of long waves. 
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Appendix A: base flow coefficients 

Here we detail the coefficients, (j)i{y] xq, Re) and 
Xiiu', Xq, Re), of tfie asymptotic expansion representing 
the intermediate and far base flow. This approxima- 
tion is homogeneous in the x and z directions, and pa- 
rameterized through the downstream station xq and the 
Reynolds number Re. 

Zero order, i=0 

(/-o = Co (Al) 
Xo = (A2) 

with Co = I. 

First order, i=l 

(jji ^ -^C'ie-"'=^'/('^^«) (A3) 
XI = (A4) 

with Ci = 1. 



with Cg = -2.26605 + 0.15752 • Re - 0.00265 • Re^ + 
0.00001 • Re^. 

Coefficient A is related to the drag coefficient Cd 
(A = |;(i?e/7r)^/^C£)(i?e), iFi is the confluent hyperge- 
ometric function, Hr„_i(?7) = H„_i (ii?e^/^77), where ff„ 
are Hermite polynomials, and 

Re 2 

Fniv) = / ^V' , G»(^)d7?; (A9) 
Gniv) - A~'' J AUij)m-n-i{v)dv ■ (AlO) 
where 77 = y/ ^/xq. 



Second order, i—2 

+ -^V^erf(-J— y)] (A5) 
2 va;o 2 V a;o 

X2 = -4^e-«^^^/(^-) (A6) 

with C2 = -2.75833 + 0.21237 • Re ~ 0.00353 • Re'^ + 
0.00002 •i?e3. 



Third order, i=3 



</,3 = A3e-^^^'/('*"«)(2-i?e^) 



x[-C3 - ReF3{xa,y)] 
2 2 ^/xo^ 



(A7) 



X / [c-^<'/(4^«)iFi(- 



1 1 i?eC^ 







to 



1 ^e-«e.V(4.o)^F^(_l,l.:^)] 



2 y/x^ 

_}__y_^-B^ey^/{2xo) 

2 ^/i^ 



2 V Xo 



(A8) 



Appendix B: ordinary differential operators of the 
initial-value problem 



In this Appendix we list the ordinary differential opera- 
tors of the full linear system and, up to order 0{k), of the 
multiscale system. The coefficients inside these operators 
are the quantities computed in xq that are associated to 
the spatial derivatives of the vorticity and velocity of the 
base flow (U=(C/(y; Xq, i?e), xo, i?e))). In particu- 
dV dU 

lar, fiz = (— 7^)\x=xn is the mean vorticity in the 

ox oy 

spanwise direction, and the coefficients a, 6, c, d, e are 
equal to: 



dx 

d~ — 
dx 



dx"^ 



dxdy 



dV 

dx 



i{kcos{4>) + iaAU - V— + 
dy 



Re^dy^ 



2ikcos[(j))ai\ + 



i(kcos((j)) + iai) d 

u + tt; —. 7-T 



k"^ + 2ikcos{(j))ai — dy^ 
iksin{(j)) 



^ks^ni^f^ 
dy 



92 



k^ + 2ikcos{(j))ai — of dy^ 



Order 0(1): 



Full linear problem: 



G 



i(kcos{(f>) + iai)U 



dy 



1 

Re dy' 



(Bl) 



H= - 



+ 



K- 



i{kcos{(j)) + iai) ^ d 
fc^ + 2ikcos((/))ai — q;| dy 

■n ■ ^^^t, 

t[kcos((p) + la,;)- 



dy 



fc2cos2(0) + 2ikcos{(t>)ai 



k'^ si-n? {(j)) 



fc^ + 2ikcos((j))ai 



9 ,2 ■ 2^^,9V 
X a— + fc'^szTi^ ((/))— + 
dy dy 

+ {k'^cos'^{(l)) + 2ikcos{(j))ai - a- )d + 



^(0) + 2ikcos{(j))c 



fc^ + 2ikcos{(j))ai 
k'^sin'^{4)) 



-af d^ 

2 



- a_ 
dV 92 



k'^ + 2ikcos{(f>)ai — af dy dy"^ 
d 

— i{kcos{(j)) + iaAe— — h 
dy 

i{kcos{(j)) + iai) d^ 
fc2 + 2ikcos((j))ai — af dy^ ' 
ksin((j)) 



fe2 -f 2ikcos{(j))ai 



b — iksin((j))e + 



^{kcos((j)) + iai)ksin{4>) 
fc2 4. 2ikcos((j))ai — af 
{kcos{(l)) + iai)ksin{(j}) ^ d 
k^ + 2ikcos{(j))ai — af dy 
{kcos{(j)) + iai)ksin{(j)) dV d 
k^ + 2ikcos{(l))ai — af dy dy 



+ 



ksin{(j)) 



d^ 



k^ + 2ikcos{4>)ai — af dy^ 



(B2) 



(B3) 



Go — 
Ho = 

Lo = 



d 



dy ai dy 



- d 



a,U -V 
1 



dy 
d 



a,^ ^ dy^ dy^'' 



1 



Re dy 



ai^ dy' 



Order 0(k): 



Gi 



-icos{(l))U - V 



d 



1 a' 

+ — 2 

Re} dydY 



dY 

2icos{(j))ai], 



9 1.5 
dY ai dY 

- —cos{(j))b- icos[(j))—— 

a^ dy dy 



d^ 



dydY 



d + 2iaiCos{(j))d + 



3 d^ d_ 



Q3 



-^cos(4})e—-^ — icos((f))e 
t oy^ 



sin{(f))a ^sin{(j))b 



d_ 

dy' 



Oil 

isin{(j))e{l 



af dy"^ 
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d 

Li = -icos(0)[/ - y— + 



Ml = :^sin[(l))e—-^ . 

ay af oy^ 
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